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Magnetic fields in planetary nebulae and post-AGB nebulae 
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ABSTRACT 



X 



Magnetic fields are an important but largely unknown ingredient of planetary nebulae. 
They have been detected in oxygen-rich AGB and post-AGB stars, and may play a role in the 
shaping of their nebulae. Here we present SCUBA sub-millimeter polarimetric observations 
of four bipolar planetary nebulae and post-AGB stars, including two oxygen-rich and two 
carbon-rich nebulae, to determine the geometry of the magnetic field by dust alignment. Three 
of the four sources (NGC 7027, NGC 6537 and NGC 6302) present a well-defined toroidal 
magnetic field oriented along their equatorial torus or disk. NGC 6302 may also show field 
lines along the bipolar outflow. CRL 2688 shows a complex field structure, where part of the 
field aligns with the torus, whilst an other part approximately aligns with the polar outflow. 
It also presents marked asymmetries in its magnetic structure. NGC 7027 shows evidence for 
a disorganized field in the south-west corner, where the SCUBA shows an indication for an 
outflow. The findings show a clear correlation between field orientation and nebular structure. 
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1 INTRODUCTION 



Whether magnetic fields play a role in the shaping of planetary 
nebulae (PNe) is an open question. Most of the po st-AGB nebulae 
appea r elliptical, bipolar or even multi-polar dBalick & Frankl 
l2002h . These morphologies are amplified by the interaction 
between a slow AGB wind with a faster post-AGB wind. However, 
this amplification still requires an initial asymmetry in the slow 
wind. This initial shaping has been mainly attributed to two 
possible phenomena: binarity and magnetic fields. 

In the binary model, a close companion affe cts the mass-losing 
AGB star via common envelope evolution jPe Marco & Mod 
2006), mass transfer and/or tidal forces, and may lead to the 
formation of an accretion disk around the companion star. The 
binary orbit provides a source of angular momentum which 
get carried by the wind, and leads to an equatorial disk. The 
angular momentum loss by the stars may lead to a merger. The 
model of shaping due to this binary interaction is quite popular 
jBond & Murdini |2000|) . dCiardullo etafl 120051) . partly because 
of the presumed impossibility for a single star to supply the 
energy nece ssary to crea te a magnetic field strong enough for 
its shaping dSoked 120051) . Nevertheless, there is still a lack of 
observational evidence for the occurrence o f close binary sys tems 
during the AGB phase. Neither companions (Riera et al. 2003) nor 
high orbital velocities of the AGB stars (Bar nbaum et alj 1 19950 
are detected in a sufficie nt amount to es tablish the role of the 
binarity as predominant jMatt et al] |2000). Evidence for binary 



interactions may be found in 25-50 per cent of planetary nebulae 
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(Zijlstra 2006), althou gh values of 50-100 per cent have a lso been 
suggested dDe Marco et al.ll2004lDe Marco & M oe 2006). 

On the other hand, the magnetic field may act as a "squeezer" 
around the central star of the PN and thereby give the dust its di- 
rection (in the sense of th e outflow). Magnetic fiel ds have been de- 
tected arou nd AGB stars jVlemmings et al]|2006t) and a few post- 
AGB stars jBains et alj|2004) using radio observations of masers 
(H2O, SiO, OH). Such observations measure a local value of the 
strength of field, within masing high-density clumps, which may 
differ from the global field. The origin of the magnetic field is 
unknown: a dynamo effect resulting from an interaction between 
a slow rotating envelope an d a fast rotating core has been pro- 
posed (Blackman et al. 2001). As magnetic fields are now known 
to be present in the AGB and the post-AGB phase, their importance 
should not be ruled out prematurely. 

iGreavesI (2002) found evidence for dust alignment in two 
carbon-rich objects, NGC 7027 and CRL2688 which are respec- 
tively a young and compact bipolar PN and a strongly bipolar proto- 
planetary nebula (PPN). This was based on polarimetric Scuba 850- 
micron observations. The data suggests the presence of toroidal 
collimated magnetic fields, as would be required for the shaping. 
But the presence of such a field was not conclusively proved, be- 
cause of the very few detected vectors and limited spatial resolu- 
tion. We present here new 850-microns and the first 450-microns 
polarimetry of post-AGB stars, which allows for better resolution 
and reduces angular smearing. In addition to the objects observed 



2 L. Sabin et al. 



bv lGreavesI J2002t) . we also observed NGC 6537 and NGC 6302. 
The sample contains two carbon-rich and two oxygen-rich nebulae. 
We show the presence of well-aligned toroidal fields in three of the 
nebulae. The fourth object shows indications for both a toroidal 
field and one aligned with the polar outflow. We conclude that 
toroidal fields may be common in bipolar PNe, and could play a 
role in the shaping of the nebulae. 



2 OBSERVATIONS 

The polarimetric data have been obtained May 10 th 2005, with the 
polarimeter on the Sub-millimeter Common-User Bolometer Array 
(SCUBA), a t the JCMT. The inst rument (now decommissioned) is 
described in lHolfandetai1dl999h . The JCMT beam size is 15 arc- 
sec at 850/im and 8 arcsec at 450 /im. 

SCUBA contains two arrays of bolometric detectors, covering 
a field of view of 2.3' in diameter. The 850-/im array has 37 individ- 
ual detectors, and the 450-/im array has 91 detectors. The two ar- 
rays are used simultaneously. The spatial resolution of the detector 
is 7.5" at 450 /jm, and 14" at 850 /im. The gaps between the detec- 
tors are covered by moving the telescope in sub-pixel steps, in the 
so-called jiggle-map mode. The step siz e needs to be optimiz ed for 
the wavelength used. The polarimeter dGreaves et alj|2003h mea- 
sures the linear polarization by rotating the half-waveplate in 16 
steps. We used this in combination with jiggle map mode. A chop- 
throw of 45 arcsec was used. 

Good photometric images cannot be obtained at both wave- 
lengths simultaneously. When the jiggle pattern is optimized one 
wave-band, the simultaneous image obtained in the other band is 
under-sampled. The size of the pixels in the final image is set 
during the data reduction with the polarimetric package of Star- 
link. For a jiggle pattern corresponding to a 450/im measurement, 
the pixel spacing becomes 3 arcsec and for 850/xm it is 6 arc- 
sec. The reduction of i nfrared polarimetry is discussed in detail by 
Iffildebrand eTaO j2OO0l) . 

The instrumental polarization (IP) of each detector has to 
be removed in order to have a correct polarization calibration. A 
new IP calibration, provided on site, was used: this contains a re- 
measurement of the central detector. For the other detectors, an 
older IP calibration was used. The IP should not vary too much 
for detectors close to the center of the array. Some tests done to 
see the importance of the instrumental polarization (by changing 
its value by the size of the error for example ) and if it could affect 
the polarization vectors length, showed that in our case, the IP was 
small (at850^m: ~1.20%±0.25%andat450/xm: ~3.25%±0.25% 
) and didn't play any role in modifying the vectors. The resulting 
polarimetric images were checked for different pixel binning, from 
lxl to lOx 10, and we didn't find any significant changes. 

The linear polarization is measured as a percentage polariza- 
tion, and a direction. The polarization is typically caused by the 
alignment of spinning dust gr ains, with their long axis perpendicu- 
lar to the local magnetic field dGreaves et al. 1999). Thus, the mea- 
sured angle of polarization is 90 degrees rotated with respect to 
the magnetic field. The degree of polarization does not give direct 
information on the strength of the magnetic field. 

We observed four targets: NGC 6537, CRL 2688, NGC 6302 
and NGC 7027. They were observed in jiggle map observing mode 
at both wave-bands, 450^im and 850^im. For each object, the mean 
direction and angle of polarization are listed in Table 1 . 

The figures below show the direction of the magnetic field, i.e. 
the vectors are perpendicular to the direction of the grain alignment. 



3 THE RESULTS 
3.1 NGC 6537 

NGC 6537 is also known as the Red Spider nebula. It is a bipolar 
planetary nebula with a very hot central star (1.5-2.5 x 10 J K). The 
nebula suffers extinction by circumstellar dust. This extinction is 
locali zed mainly in a comp act circumstellar shell. An extinction 
map dMatsuura et alj [20051) reveals a compact dust shell with a 
roughly spherical inner radius of about 3 arcsec, with a minimum 
towards the central star dMatsuura et ai] |2005l: ICuesta et aljfl995h . 
The pola r outflows (traced by hi gh velocity winds of about 
300km/s: ICorradi & Schwarj 1 19931) extend 2 arcminutes along 
the NE-SW direction. The Scuba 850-^tm continuum map (Fig. 
1) shows elongation perpendicular to the outflow direction, with 
a major-axis diameter of approximately 20 arcsec. The extinction 
map is now seen to represent the inner edge of a more extended, 
and possibly toroidal, structure. The equatorial plane may be 
oriented a little closer to the EW direction, based on the Scuba 
map. 

The best polarimetric results obtained with SCUBA are ob- 
tained at 850/xm. The consistent orientation of the polarization vec- 
tors shows that the magnetic field (hereafter B) has a dominant 
direction along the equatorial plane, approximately perpendicular 
to the outflow direction. The dust alignment is therefore directed in 
the same sense as the outflow. The length of the 18 different vectors 
does not show large variations, with a degree of polarization vary- 
ing from 8 to 14%, suggesting a consistent magnetic field. More- 
over the absence of smaller polarization vectors toward the center 
of the nebula indi cates that there is no change in geometry of B 
towards the core dGreaves! I2OO2I) . This supports a location some 
distance from the star (i.e. in a detached shell), since otherwise av- 
eraging of vectors in different directions within the JCMT beam 
would reduce the detected net polarization at the central position 
(beam depolarization). 

These observations indicate the presence of a consistent 
toroidal magnetic field, located along the equatorial plane of NGC 
6537, in a circumstellar torus. Compared to the size of the outflow 
lobes, the field is located relatively close to the star. The presence 
of a B-field had already been su spected, based on the occ urrence 
of filaments near the central star (Huggins & Manlev 2005). 

The extinction map obtained bv lMatsuura et alJd2005h (Fig. 1- 
Bottom panel) shows an asymmetry, in that the extinction is higher 
on the western side (^2 mag versus £1.6 mag towards the east). 
This asymmetry is also seen in the Scuba data, with a larger ex- 
tension and more polarization vectors on this side. This is a further 
indication that the magnetic field is located within the detached dust 
shell. 

There is no strong indication for magnetic fields along the spi- 
der lobes. There is a slight trend for the vectors to curve, but this 
is caused by only a few of the vectors and would need confirma- 
tion. We cannot state with confidence whether the lobes also carry 
a magnetic field. 



3.2 NGC 7027 

The 450/xm jiggle map of the young planetary nebula NGC 7027 
(Fig. 2) and its near environment covers a field of about 40x36 
arcsec 2 . The central star is surrounded by an ionized area of ~283 
arcsec 2 , which is in turn enclosed by a thin atomic and molecu- 
lar layer that is seen for instance in the H2 emission observed by 
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Figure 1. Scuba 850/im results on NGC 6537. North on the top and East on 
the left. The axes give the image scale in arcseconds. Top panel: the 850^m 
continuum map of NGC 6537, with contours at 1%, 2% ,5% and 10% of 
the peak. Middle panel: Magnetic field orientation. The general outflow di- 
rection of the nebula is indicated by A and the equatorial plane by B. The 
polarization vector scale on the left is set at 10%. Bottom panel: Extinction 
map presented by Matsuura et al. (2005). The highest levels of extinction 
occur at ~ 4 arcsec from the central star with E(H/3 — Ha) > 2 mag, co- 
incident with the dust emission in the 850/im map. The polarization vector 
scale is set at 10%. The magnetic field which coincides with the area of in- 
ternal extinction, is mostly aligned along the equatorial plane, indicating a 
toroidal field. The dust alignment is perpendicular to the vectors displayed. 




Figure 2. Scuba 450/ttm results on NGC 7027. North on the top and East 
on the left. The axes give the image scale in arcseconds. Top panel: The 
450^tm continuum map of NGC 7027. Contours are set from 1 to 5% and 
10% of the peak. Middle panel: Magnetic field distribution. The general 
outflow direction is indicated by A and the equatorial plane by B. The po- 
larization vector scale (showing the degree of polarization) is set at 10%. 
Bottom panel: The Scuba polarization vectors are shown on an continuum 
subtracted H2 (color-inverted) map of NGC 7027 (North on the top and East 
on the left). The field is mostly directed along the equatorial plane. The cen- 
tral region shows much reduced polarization and the field orientation differs 
on the extreme western side. 
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ICox et al.1 J20021) . A thin dark ring delineates the equatorial plane 
in optical images. The structure is better seen in the HST image 
(Fig. 2-bottom), which shows both t he surrounding mol ecular layer 
(in H2) and the inner ionized region. lLatter et al J (12000) shows that 
the NW lobe is blue-shifted (closer) and the SW lobe is red-shifted. 

The dominant direction of the magnetic field coincides with 
the equatorial plane. But this behaviour is mainly seen on the 
North-East side while the south-west part seems to be disturbed: 
the magnetic field may be "broken". The degree of polarization is 
8.9%±0.9%) towards the NE direction (or lobe) and 7.6%± 1.3%) 
towards the SW. The degree of polarisation is much reduced in the 
center of the nebu la and lacks a co herent direction here. This effect 
was also noted bv lGreavesI {2002) and may indicate that coherence 
is lost in the ionized region. 

3.3 CRL 2688 

The proto-planetary nebula CRL 2688 (or Egg Nebula) is a bright 
carbon-rich bipolar object characterized by two pairs of searchlight 
beams superposed on a reflection nebula. The origin of the light 
beams has been suggested to be a very close stellar companion 
( Sahai et al ] | 19981 ; iKastner & Sokerj|2004|) . or the presence o f dust 
layers reflecting the light of the central star dGoto et al. 2002). The 
bipolar reflection nebula shows a dark equatorial lane where a large 
amount of dust may by present. 

The 850- ^im is presented in Fig. 3-Top panels. It shows elon- 
gation in the direction of the polar outflows. The equatorial extent 
is less but at the lowest contour the torus is more extended on the 
eastern side. The extent is ~60"x40" although it is not clear pre- 
cisely where the emission ends. The bright core at 850/tm shows a 
FWHM of 16 x 14 arcsec 2 , at a position angle of 25 degrees which 
is in the same direction as the outflow. We could not measure the 
core elongation at 450/xm due to the under-sampling. The spatial 
resolution is insufficient to separate the two light beams, but the 
width of the emission suggests that the dust traces the larger lobes 
which the light beams illuminate, and that these beams themselves 
are not present in the far-infrared data. (The elongation of the core 
is in fact along the line of one of the two beams only.) 

CRL 2688 presents the highest number of polarization vectors 
in our sample. For this nebula, the under-sampled 450/im map gives 
some additional and useful indications (even if we observe a lack of 
vectors due to this under-sampling). The polarization vectors cover 
the full nebular extent (seen at both wavelengths), as they do in 
NGC 7027. The degree of polarization is not uniform and strongly 
decreases near the center. This phenomenon is more visible in the 
850/iin map, which suggests that beam depolarization may play 
a role. At 850 /im, the mean value of polarization of the region 
containing the bipolar lobes is about 3.2%, the outer region has a 
mean value of 8.8%, and the region of the dark lane shows a mean 
degree of polarization of 1.4%. 

If we draw a line passing through the nebula in a longitudinal 
direction (NNE-SSW fig. 3 bottom-left panel), we can see distinct 
behaviours of the magnetic field on either side of this line: on the 
eastern side, the orientation of B is mainly in the direction of this 
line, while on the western side, the magnetic field appears perpen- 
dicular to it. 

The 450-/im map (Fig. 3-bottom-right panel) shows the mag- 
netic field at higher resolution, confirming the bimodal distribution. 
This map is under-sampled and only some positions in the nebula 
are covered. The map gives a suggestion of a superposition of a 
toroidal field and one aligned with the polar outflows. The field 
becomes less ordered towards the tip of the outflow direction, but 



the emission here is faint and the uncertainties on the polarization 
vectors are larger. 

The complicated magnetic morphology makes it likely that the 
dynamics have shaped significant parts of the field, rather than the 
magnetic field shaping the nebula. The structure may show the su- 
perposition of two components over most of the area of the source. 
The lack of polarization in the centre suggests that the bright core 
is not polarized, or its polarization is averaged out over the JCMT 
resolution. The outermost field vectors in the longitudinal direction 
are along the outflow, and this may indicate a field carried along 
by the outflow. This is most apparent in the north when looking at 
the outer vectors direction in the figure 3-bottom-left, in the south 
the outer vectors point towards the corresponding outflow so we 
assume that the magnetic field may be carried by this outflow. 

The grain alignment, which is perpendicular to the direction of 
the magnetic field, is in the direction of the outflows on the western 
side. On the eastern side, the grain alignment is toroidal. 

High resolution molecular image s in the literature include 
IRA M CO J=2-l data jCox et al.ll2000h and HST infrared H 2 im- 
ages dSahai et alj|l998l) . The continuum-subtracted H2 image in 
Fig. 4, shows multiple jets, located both in the equatorial plane 
and towards the polar directions. The CO data show these jets to 
be the tip of flows originating much closer to the star. These jets 
all fall within the Scuba 850-/im core, indicating that indeed struc- 
ture is present on scales much smaller than the JCMT beam size. 
Thus, beam dilution and beam-depolarization of the magnetic field 
is likely. 

We note that the elongation visible in the Scuba image, on the 
eastern side, is in the sa me direction as th e largest molecular sub- 
jet E2 in Fig.4, or Al in ICoxet all (2000) (but the Scuba structure 
is much larger). The polarization vectors seem unaffected. On the 
western side, both the molecular and the Scuba images indicate a 
smaller extent of the envelope. The polar outflows as seen in the 
molecular lines do not show the light beams, but instead show rel- 
atively well-collimated lobes. The lobes are brighter in the north. 
The polar lobes in the Scuba image are consistent with this, both 
in direction and in brightness, but again are much larger than seen 
in the molecules. The molecular emission shows the wind-blown 
cavities, while the Scuba emission arises in the surrounding shells. 

3.4 NGC 6302 

This oxygen-rich planetary nebula jPottasch & Beintemalll999T) is 
a proto-typical butterfly-type nebula. It shows two extended lobes 
that result from a bipol ar outflow, with a d ark la n e in the center. 
It has been studied by IMeabum & Walshl Jl980h. IMeabum et al.l 
( l2005l) . lMatsuura et alj ( 120051) and lCasassus et alifeOOOl) . 

The Scuba observations were optimized for 450-/im. They 
show a bright core, of FWHM 14 x 12 arcsec 2 (larger than the 
beam), with the long axis in the north-south direction, with an ex- 
tension towards the south-southeast (Fig.5). We underline the fact 
that the long horizontal line in the eastern part of the nebula and 
the long vertical line in the South- West (Fig.5 bottom panel), are 
artifacts from the HST image and not Scuba polarization vectors. 
The arcs seen at 25 arcsec are likely part of the diffraction pattern 
from the very bright core. The polarimetric data obtained at 450/im 
shows only five polarization vectors. They do not line up with either 
the dark lane or the outflow. However, they are fairly well aligned 
with the ellipsoidal radio source in the center of the nebula. The 
radio nebula shows the inner ionized region, confined by the dense 
torus: the elongation is perpendicular to the torus. 

The difference in position angle of the inner torus with the 
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Figure 3. Scuba 850 and 450/xm results on CRL2688. North on the top and East on the left. The axes give the image scale in arcseconds. Top left panel: The 
850/im map of CRL 2688, with contours from 1 to 5% and 10% of the peak. Top right panel: The 850-/xm vector polarization showing the magnetic field 
orientation. The general outflow direction is indicated by A and the equatorial plane by B. The polarization vector scale (showing the degree of polarization) 
is set at 10%. Bottom panels: HST images (WFPC2, filter F606W) with the overlaid magnetic field corresponding to the 850/^m map (left) and 450/^m 
map (right). The dust continuum is elongated along the outflows. The field shows a complicated structure with components along the outflow and along the 
equatorial plane. The polarization vector scale (showing the degree of polarization) is set at 20% for both maps. 



dark lane is interpreted as a warped disk (Matsuura et al.ll2005h . 
Further from the center the outflows have different position angles, 
and eventually become east-west in the outer bipolar lobes. The 
450-ftm polarization indicates that the magnetic field may be ori- 
ented in the direction of the inner outflow. 



We have no data for magnetic fields elsewhere in the outflows 
even at 85Ctyim. 



4 DISCUSSION AND CONCLUSION 

The sub-millimeter observations of the four post-AGB objects re- 
veal new information regarding the distribution of the magnetic 
fields, the dust emission and the link between the two components. 

Extended dust emission was seen for all four nebulae. For 
CRL 2688, the large polar lobes have been detected. The orien- 
tation of its polar lobes is along only one of the two light beams. 
The equatorial emission is also extended, and shows an asymmetry 
which is correlated with the appearance of the (much more com- 
pact) equatorial jets. These jets originate close to the star and it is 
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Figure 4. Continuum-subtracted H2 (HST) image showing the jets in CRL 
2688. [North at the top and East at the left]. They are named as in Sahai et al 
(1998). The jets are all located within the compact 850-/^m source, and the 
some perturbations seen in the Scuba data can be related to the "sub-jets" 
E2or Al. 



Name 


Band (pm) 


Deg (%) 


Angle (deg) 


NGC 6537 


850 


11.2±2.2 


26.5±5.7 


NGC 6302 


450 


11.4±1.6 


32.7±4.6 


NGC 7027 


450 


8.2±0.9 


-18.5±3.7 


CRL2688 


850 


6.8±1.0 


-22.6±4.5 



Table 1. Mean values of the polarimetric parameters for the Post-AGB ob- 
jects, Deg: degree of polarization, Angle: position angle of the polarization. 
These values concern the main band studied for each nebula. 

unlikely that they are affected by structures at much larger scales. 
Instead, the jets may affect the dust emission through heating. In 
NGC 6537, an asymmetry in the dust correlates with differences 
in the extinction map, suggesting that here the Scuba map shows a 
true asymmetry in the dust distribution. 

The presence of detectable polarization shows that at least 
some of the dust grains are not spherical. 

4.1 Targets 

Polarization is detected for all four nebulae. This in itself suggests 
that magnetic fields are common for these types of objects. It is thus 
important to discuss what types of objects were selected. 

First, the four targets are all bipolar nebulae. The least pro- 
nounced morphology is shown by NGC 7027, but even here the 
molecular maps clearly show the underlying bipolar structure. We 
do not have information on the existence of magnetic fields in less 
bipolar (e.g. elliptical or spherical) nebulae. 

Second, although three of the nebulae have very hot central 
stars and ionized cores, all have extended molecular envelopes. It 
is argued above that the detected fields are located in the molecu- 
lar or atomic regions. The presence of molecular envelopes around 
ionizing stars indicates very dense nebulae, and therefore very high 
mass loss rates on the AGB. 

Third, the dense nebulae with relatively small ionized regions 
imply that their stars have evolved to high temperatures before the 




Figure 5. Scuba 450/^m results on NGC 6302. North on the top and East 
on the left. Top panel: The map shows the magnetic field orientation on the 
planetary nebula with contours from 1 to 5%, 10%, 20% and 25% of the 
peak. The polarization vector scale (showing the degree of polarization) is 
set at 20%. Note that the arcs seen 20 arcsec out are part of the beam side 
lobes of the very bright core. Bottom panel: 450/im Scuba map of NGC 
6302 overlaid on the F656N band HST image (WFPC2). North on the top 
and East on the left. The external hexagonal contours do not belong to the 
star but are representative of the bolometers. The main Eastern outflow seen 
on the submillimeter image coincides with the one observed with the HST. 
Note that the arcs seen 25 arcsec out are part of the beam side lobes and 
the two long lines (horizontal and vertical) on both sides of the nebula, are 
artifacts from the HST image and not polarization vectors. 

nebulae have had time to expand significantly. Such fast evolu- 
tion is characteristic for high-mass central stars. The cooler object 
(CRL 2688) may have a lower mass, but the fact that it has become 
a carbon star suggests it still has a relatively higher initial mass 
{Mi <<5 3Mq), as lower-mass stars experience insufficient dredge- 
upo to become carbon rich. For NGC 6302, a progenitor mass of 
4-5 Mq has been suggested, base d on the mass of the circumstel- 
lar envelope l lMatsuura et alj|2005t) . The precise initial masses are 
not well determined, but high mass progenitors appear likely for all 
four objects. 

Thus, the result of the survey can be interpreted that magnetic 
fields tend to be present for high mass progenitors evolving into 
bipolar nebulae, and that the fields are detectable while the nebulae 
still have molecular envelopes. 
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4.2 Field location and orientation 

Overall, we see fields aligned with the polar directions and/or 
toroidal fields in the equatorial plane. Assuming that the dust in the 
polar flows originates from the dust reservoir in the disk or torus, 
this suggests that the field also originates there. The alignment with 
the polar flow is interpreted that the gas carries the field with it. 

NGC 6302 is the only source in our sample without evidence 
for a toroidal field. This nebula has a complicated morphology with 
a warped disk and the nebula is oriented at 45 degrees with respect 
to the equatorial dust lane. It is therefore difficult to ascertain the 
alignment of the field with any particular morphological compo- 
nent. However, alignment with the radio core seems most likely. 

The other nebulae show evidence for toroidal field compo- 
nents. We therefore suggest that the initial configuration of the field 
is toroidal, and is located in the equatorial torus. 

The polarization in all cases appears to be seen in the neu- 
tral/molecular regions. Towards the central, ionized regions, the 
percentage of linear polarization is much reduced. This may in part 
be because of the beam averaging over regions of different polar- 
izations, but the ionized regions may also be expected to lose a 
directional field (non-constancy of B). The reason for the absence 
of a detectable magnetic field in the ionized region can also be ex- 
plained by the lack of dust (so B is not carried). For NGC 6302, the 
field is aligned with the ionized core but is likely located around 
rather than within the core. 

Polarization will not be detected if the field is oriented along 
the line of sight. This is expected for a toroidal fields near the tan- 
gential points of the torus. This may be present especially in CRL 
2688, where two components are found at different locations. 

4.3 Chemistry 

The sample includes two oxygen-rich nebulae (NGC 6537 and 
NGC 6302, both PNe with hot central stars) and two carbon-rich 
nebulae (NGC 7027 and CRL 2688: the latter is the cooler post- 
AGB nebula in our sample). This sample is too small to draw gen- 
eral conclusions. However, we point out the following points of 
interest, for future study. 

• The O-rich PNe both present collimated magnetic fields con- 
centrated near the central star. In both cases, we have an organized 
magnetic field in the central region, resistant to distortion effects, 
but do not detect more distant fields. 

• The two carbon-rich nebulae show extended fields covering 
their entire nebula. In both nebulae, the magnetic field is more dis- 
organized than in the two O-rich stars. 

The question whether there is a chemistry-related origin to 
these behaviours is valid, but cannot be answered with the cur- 
rent sample. Dust chemistry could play a role. First, carbon grains 
tend to be smaller than silicate grains. Larger grains may be more 
spherical, which would give less (or no) polari zation. Both oxygen- 
rich nebulae show evidence for crystalline ice jMolster et alj2002l) 
(strongest in NGC 6302), indicating coatings on the grains. Second, 
amorphous carbon grains may be intrinsically non-spherical, build 
up from graphitic carbon sheets. 

For the O-rich nebulae, it is possible that only the hotter 
dust, above the ice evaporation temperature, leads to non-spherical 
grains and a polarized signal. Thus, a difference in polarization 
does not necessary imply a difference in magnetic field structure. 

Both O-rich nebulae show evidence for weak PAH emission in 
their central regions. The origins of the PAHs is not clear, but they 



are seen only in irradiated regions. The continuum radiation from 
the dust component containing the PAHs is only seen at short wave- 
lengths (warm dust) and contributes very l ittle to the sub-mm flux 
dKemper et alj|2002l ; iMatsuura et alj|2005f) . Thus, the polarization 
signal is seen in the oxygen-rich grains, not the PAHs. 

Iron needles could be considered as carriers of the polariza- 
tion: in O-rich stars, some or most of t he iron is incorporated in the 
silicate grains (olivines and pyroxenes: iFerrarotti & Gail 200j]), but 
for a C/O ratio close to an d above unity, solid ir on and FeSi be- 
come more important (e.g., Ferrarotti & Gail 2002). Metallic, non- 
spherical iron grains may contribute to th e NIR opacity in high 
mass-loss OH/IR stars ( Ke mper et alj|2002l) . 

4.4 Comparison less evolved/ more evolved nebulae 

The extent of the nebulae shows that the objects differ in age. The 
nebulae of CRL 2688 (Post-AGB) and NGC 7027 (young PN) are 
younger than NGC 6537 (PN) and NGC 6302 (PN). (The lobes of 
the last object extend to over 2 parsec across.) 

Our submillimeter polarimetric data show that the magnetic 
field is better organized for the older nebulae. Assuming an orig- 
inally toroidal field, later carried by the bipolar flows, the best 
organized field is expected if one component dominates. The 
bipolar flows are mos t massive early in the post-AGB evolution 
flBuiarrabal et alJuOOlh . As they diminish, the remaining field may 
become less confused. However, evolving dust characteristics may 
also play a role: dust grains become larger in older disks, and may 
show less polarization as a result. It is not possible with the current 
sample to separate evolutionary and chemical effects, but both are 
expected to occur. 

4.5 Field strength and origin 

The detected polarization contains information on the direction 
of the field, but not its strength. A method to obtain the field 
strength out of fluctua t ions in polarization was suggested by 
Chandrasek har & Fermil dl953l) . Applying this method to our data 
leads to field strengths of the order of mG. This equation assumes 
that the magnetic field dispersion is caused by Alfven waves. The 
applicability to our sources may be in doubt, as the dispersion may 
be dominated streaming motions. 

In neutral-dominated media, Alfven waves contribute if the 
collision times between ions and neutrals are shorter than the period 
of the Alfven wave. This is the case for waves longer than 




where n is the particle den sity per cm , and /; is the fractional ion- 
ization (Hildebrand 1996). For our sources, /, may be high except 
in the densest regions of the tori. Assuming a typical source size of 
0.05 pc, and a field strength of a mG, Alfven waves may contribute 
if fi > 10~ 4 . This is a plausible value for the outflows. 

Previous observations of OH Zeeman splitting has shown field 
strengths in OH/IR stars and post-AG B stars of a few mG a t 
3 x 10 15 - 2 x 10 16 cm from the stars dBains et al.|[2004 l2003h . 
Much larger values have been inferred from wa ter masers: several 
Gauss at ~ 2 x 10 14 cm jVlemmings et al.l2005l) . The difference is 
consistent with a dipole field, which gives an r~ 3 dependence, and 
less consistent with a solar-type field r~ 2 . Water masers trace high 
density clumps and the measured field strengths may not be fully 
representative of the surrounding areas. 
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Figure 6. . Location of the Scuba dust emission, compared to the HST images (contours). We present respectively NGC 6537, NGC 7027 and CRL 2688. 
NGC 6302 is shown in Fig. 5 (bottom panel). 



Assuming we look at typical distances of 5 arcsec at 1 kpc, 
our measurements are at ~ 5 x 10 16 cm from the stars. Compared 
to the OH and water masers, one would expect a field of lmG or 
less at this distance. This suggest that the direction changes in the 
field are due to streaming motions, rather than Alfven waves. 

The grains become aligned with the magnetic field with the 
Davis-Greenste in mechanism. The time scale for this to occur is 
dKruegefeOOi Ch. 11) 

t re \ oc B~ 2 

and, for fields of the order of a few mG, is typically 10 6 yr. The 
age of the nebulae is ~ 10 4 yr. For this time scale, the alignment 
of the dust with the magnetic field requires fields in excess of those 
seen from OH maser emission. The alignment is therefore likely to 
originate close to the central star, at r < 10 15 cm, and is maintained 
while the nebula expands. 



4.6 Evolution 

For OH /IR stars, evidenc e for d ipole fields has been presented 
by, e.g.. lVlemmingsetalJ EM). For more evolved stars, in the 
early post-AGB evolution, the presence of both toroidal (IRAS 
20406+2953) an d poloidal (OH17 . 7-2.0 ) fields are inferred from 
OH observations (Ba ins et alj2004l . l2003h . The poloidal fields arise 
from a stretched dipole. For the even later evolutionary stage stud- 
ied here, a strong toroidal field component is found, combined with 
a poloidal field. 

A toroidal field can form out of a dipole field by rotation. This 
already makes it likely that the formation of the torus and the for- 
mation of the toroidal field are related. A binary companion can 
be the source of the required angular momentum. Rotation of the 
star itself is less likely, as its angular momentum can be expected 
to be lost early on in the mass loss. A binary companion, on the 
other hand, can deposit its angular momentum during the time of 
the peak superwind, at the end of the AGB mass loss. 

The fields detected via OH and H2O masers have been claimed 
to be strong enough to dominate the dy n amics of the nebula 
jBains et alTl2004l ; Iviemmings et al J 120051) . ISokeJ d2006h argues 
that fields of this magnitude cannot be generated by the star itself, 
and should be attributed to a companion. Based on this, the basic 
shaping mechanism of the nebula is found in binarity. However, 
once the strong fields have been generated, they would be a signif- 



icant factor in the further evolution of the nebula towards the PN 
phase. 



4.7 Summary 

We present the discovery of magnetic fields in four bipolar post- 
AGB stars (NGC 6537, NGC 7 027, NGC 6302 and CRL 2688). 
This confirms the earlier work o f G reave si J2002I) for two of these. 
The fields are mapped at high resolution (for sub-mm), using either 
450 or 850/im. The sub-mm emission traces the extended emission. 
In CRL 2688, we find evidence for the polar lobes which contain 
the well-known search beams. 

All objects show polarization indicative of grain alignment 
by magnetic fields. Toroidal fields are found for three objects, and 
poloidal fields for two (CRL 2688 shows evidence for both). The 
alignment of the field with nebula is least certain for NGC 6302, 
where the nebula shows a multipolar structure. Our results sug- 
gests that magnetic fields are common in these types of targets: 
bipolar nebulae with intermediate-mass progenitors. The fields are 
long lived as they are observed over different evolutionary stages. 

The data also show evidence for elongated grains. The po- 
larization is stronger and more extended for the carbon-rich nebu- 
lae. This may show that amorphous carbon grains are intrinsically 
non-spherical, e.g. the sheet-like structure of graphite. The oxygen- 
rich nebulae show polarization only in their central regions. This is 
interpreted in that the dust grains in the torii are larger and more 
spherical. 

The poloidal fields in the polar flows suggest that here the field 
is carried along by the flow and that these flows are not magneti- 
cally confined. The toroidal components in the torus may be more 
important, for the dynamics. As the fields in AGB stars are dipole- 
like, the toroidal field structure appears to be a later phase of evo- 
lution. We suggest that this transition occurs when the equatorial 
field is wound up through the interaction with a companion. Fol- 
lowing ISokeJ J20061) . we suggest that the original shaping agent 
is a binary companion. Once the field has been wound up, it may 
however become dynamically important for the subsequent nebular 
evolution. 
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